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Abstract Observations from multiple spacecraft show that there are energy spectral 
“breaks” at 1-lOMeV in some large CME-driven shocks. However, numerical models 
can hardly simulate this property due to high computational expense. The present paper 
focuses on analyzing these energy spectral “breaks” by Monte Carlo particle simulations 
of an isolated CME-driven shock. Taking the Dec 14 2006 CME-driven shock as an ex¬ 
ample, we investigate the formation of this energy spectral property. Eor this purpose, we 
apply different values for the scattering time in our isolated shock model to obtain the 
highest energy “tails”, which can potentially exceed the “break” energy range. However, 
we have not found the highest energy “tails” beyond the “break” energy range, but instead 
find that the highest energy “tails” reach saturation near the range of energy at 5MeV. So, 
we believe that there exists an energy spectral “cut off’ in an isolated shock. If there is 
no interaction with another shock, there would not be formation of the energy spectral 
“break” property. 

Key words: acceleration of particles — shock waves — Sun: coronal mass ejections 
(CMEs) — solar wind — methods: numerical 


1 INTRODUCTION 

Strong astrophysic al shocks are often associated with superther mal particle emission and with magnetic 
field amplification (iBvkov et al.ll2013[lVladimirov et alll2006h . This phenomenon suggests that shocks 
are regions where particles are efficiently accelerated, and this large group o f energetic particles is 
responsible for the excita tion of magnetic turbulence via plasma instabilities (iBelll ll978l : lBell et all 
l2013l : ljokipii. J. R'll2013h . These magnetic fields which diffuse cosmic rays in the vicinity of the shock 
are required to be much higher than the averaged magnetic field in the interstellar medium. 

The theoretical model includes the determination of the particle injection energy from the ther¬ 
mal particle distribution into the non-thermal particle distribution, the maximum energy of particles 
accelerated at the shock, energetic particle spectra at all spatial and temporal locations, and the dynam- 
ical distribution o f particles that escape upstream and downstream from the evolving shock complex 
dZank et al.Ll2000h . Monte Carlo simulation results indicate that solar ejec ta transfer energy i nto the non- 
thermal particles in an interplanetary shock with an efficiency of ^ 10% dWang et al.Ll2013l) . Studies of 
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the dependence of this efficiency on the angle between shock normal an d the magnetic field direction 
(0R/v )can have implications for ground-level enhancement (GLE) events (iLi et al. J2010l:[^odin et all 
l2013h . Estimation of the maximum particle energy by CME-driven shocks is becoming more and more 
vital for forecasting space weather. Since particles accelerated at the shock escape rather easily from the 
acceleration site, they can be detected well before the arrival of the shock. This, of course, has immedi¬ 
ate and interesting implications for space weather monitoring and prediction systems, but it does imply 
too that the study of the ion acceleration mechanism is complicated by the subsequent interplanetary 
propagation of the energetic particles. 

Eor the past several decades, there has been much literature focusing on all aspects of the diffu¬ 
sive shock acceleration (DSA). In the past, cosmic ray (CR) spectra, acceleration efficie ncy, and am 


plific ation of the magnetic field have been calculated regularly via a two-fluid appr oach dO’ C. Drury 


1983 

). More recently, those have been simulated via a particle Monte Carlo method dVladimirov et aL, 

2006 

; Wang & Yanl 201 ll Ellison et al.l 

2004 1990l Niemiec et al.l 20041). or 

na hybrid method 

(Caorioli et al.. 20141 Gamate et al.. 20121 

Giacalone et all 1993; Guo et al., 20 13HWinskell 19851), and 

via full particle-in-cell (PIC) method (Amano et al.l 2007;iRiauelme et al.. 20 

11 ). In addition, the CR’s 

transport equations have also been solved bv a numerical method dKang et al., 

2002 

Zirakashvili et al.l 

2 OIOI) and an analytical method dLiu et al.ll2004tlCaDrioli et al.Ll2010l Maikov et al.l 

1996). These meth- 


ods are all able to provide consistent results for the dynamics of the shock including the CR’s back- 
reaction. However, unlike the analytical method, the particle method a nd the numerical MH D method 
have not yet been able to simulate the energy spectral “break” property (iMalkov et al.ll20f^ . Since the 
“break” of the energy spectrum would be associated with the particle leakage mechanism, Maikov has 
presented a new combined diffusion coefficient to describe particle acceleration and escape in different 
regions. It accounts for a high turbulent magnetic field in the vicinity of the shock site (particle acceler¬ 
ation) and for faded turbulence of the magnetic field far away from the shock front (particle escape). 

Although it is widely accepted that the most efficient acceleration of solar energetic particles(SEP) 
would happen in CME-driven shocks, the underlying acceleration mechanism in the shock environ¬ 
ment still remains uncertain. In particular, it is not clear how the extensive maximum particle energy 
can be produced and why the energy spectral shape can be broken (i.e., wh y an abrupt change i n 
slope of the energy spectrum can occur ) in some large CME-driven shocks (iMewaldt et al.L l2008h . 
In the past solar cycle 23, there are several observed events exhibiting proton energy spectral “breaks”. 
These events occurred on Nov-6-97, Apr-15-01, Jan-20-05, Sep-7-05, Dec-05-06, and Dec-14-06, re¬ 
spectively. In addition, there are hard X-ray and 7 -ray energy spectra from Reuven Ramaty High Energy 
Solar Spectroscopic Imager (RHESSI) occurring on 2002 July 23. This event shows a double-power-law 
spectrum with a “break” at ~30keV in X-ray and a high energy “cut-off’ tail at ^5MeV in 7 -ray. The 
X-ray spectrum indicates that the substantial elect ron acceleration achieved to tens of keV. The 7 -ray 
line shows that ions are accelerated to tens of MeV dLin et al.Ll2003h . There are also some debates about 
broken lower energy spectrum in X-ray, which is far different from an ad hoc assumption of hot ther¬ 
mal plasma presenting as the highest low-energy cutoff (~20keV). Actually, there are a lot of analyses 
of the hardening spectra occurring in t he energy range varying from 20keV to a few MeV (iGan et al.L 
l200lUKonger^l2013t lH uangj l2009h . More recent years, an ex tensive so l ar ene rgetic particle (SEP) 
event was detected by STEREO A on July 23 2012 near 1 AU. iLiu et al.l (12014l) suggest that the in¬ 
transit interaction between two closely launched CMEs resulted in the extreme enhancement of the SEP. 
These results provide a new view crucial to space weather and solar physics as to h ow an extreme space 
weather event can be produced from an interaction between solar energetic ejecta dGonalswamv et akl 
l2005UCheng et al.Ll2013tlWang & .liLl2()l3l:ISu et al.Ll2()l3l:ISchneideRll993h . 


The parallel shocks show an effective amplification of the initial magnetic field due to the cu rrent 
of energetic ions that propagate anisotropically into the upstream flow. ICanrioli et all d20 1 3l I 20 1 4h have 
used 2D and 3D hybrid simulations with large computational boxes to reveal the formation of upstream 
filaments and cavities, which eventually trigger the Richtmeyer-Meshkov instability at the shock, and 
lead to further turbulent amplification of magnetic fields in the downstream region. The ty pical acceler¬ 
ation time, up to energy B in a shock with its velocity Vsh, is of order Tacc ~ D{E)/vlf^ dO’ C. Drurvi 
Il983h . where D{E) is the diffusion coefficient. The acceleration characteristic time scale would be 
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roughly equivalent to the time scale of the ejecta-dominated stage: when the shock is formed, the shock 
velocity drops quite rapidly, and so does amplification of the magnetic field dGuo et al.Ll2012h . In a cer¬ 
tain acceleration time scale of the shock system, the maximum particle energy is decided. However, in 
the particle simulation system, the particle’s free escaped boundary (FEB) size should be considered, 
that means the highest energy particle would be escaped from the FEB. If the present simulation will 
focus on production of the maximum particle energy, the highest energy spectral “tail” should be re¬ 
served. To obtain the maximum particle energy, we can either add the FEB size or decrease the value 
for the scattering time. Due to the expansion of the FEB size, the shock system will be brought the extra 
computational burden, so we can change the value for the scattering time. In the amplified magnetic 
field with an order of magnitude SB/Bq ^ 1, the scattering time is an important factor to determine 
the acceleration efficiency in the resonant diffusion condition, and thereby determine if the maximum 
particle energy can be achieved. In this isolated shock model, we can investigate the maximum particle 
energy by changing the value for the scattering time. 

Because we are not sure if an isolated CME-driven shock can accelerate energetic particles beyond 
Ebreafc and even up to GeV, we take an isolated shock as an example to investigate the maximum particle 
energy and energy spectral “break” by using different values for the scattering time within resonant 
diffusive scenarios. According to the DSA theory, acceleration efficiency is significantly enhanced once 
the mean free path for pitch-angle scattering is approximately equal to th e particle gyroradius (i.e. A w 
tl^E) oc E/B), and the diffusion coefficient reads Db{E) fa vrL{E) (lLagage& Cesarsk"^Il983h . If 
the Bohm diffusion condition is satisfied in the shock system and a typical interplanetary magnetic field 
with an order of a few mG, one can estimate that the maximum particle energy would be ^ 1 MeV, 
which is not enough to explain the energy spectral “break” at 1-lOMeV in observations! Ellison et al.L 
Il990h . Therefore, we hope to extensively calculate the maximum particle energy Emax using different 
values for the scattering time within an isolated shock model. If we can obtain the Emax > E^reak, 
it would imply that there is unnecessary to use multiple shocks model to explain the energy spectral 
“break” property. If we obtain the Emax < E^reak, then we should examine that whether there would 
exist an energy spectral “break” at E^reak and whether we should need a multiple shocks model indeed. 

In present paper, we do simulations to further investigate the maximum particle energy in an isolated 
CME shock by using different values for the scattering time. In section 2, we introduce the dynamical 
Monte Carlo simulation method simply. In section 3, the simulated results and analysis are presented. 
In the end, we give a summary and some conclusions. 


2 METHOD 


Many deviations of DSA arisen from the nonlinear effects of the shock, such as the modification 
of the shock structure, magnetic field obliquity, time- dependence, magnetic field amplification, and 
etc. Those have been calculated by a two-fluid model (Driirv & V61kL[l98lh and an analytical model 


dCanrioli et al.L l2010t iMalkov et all [1996Arnuto_et^ 


2006) and partic l e models including hybrid . 


particle-in-cell,and Monte-carlo method dGargate-gtall |2012 ; GiacalonS. 2004]; lAmano et al.L l200^ 
iRiauelme et all |201 ll : IVladimrrov et al.L 12006b Ellison et al. , 2004 ; Wang & Yanl 201 ih These mod¬ 
els return consistent results well, also can provide the results of the dynamics of the shock including the 
CR’s back-reactions. In general, there are two aspects of the deviation of DSA: one aspect is about the 
transfer issue depending on the macro factors of the shock in Mach number, magnetic field obliquity, 
and time-dependence etc; another aspect is about the acceleration issue depending on micro factors of 
the shock in diffusive coefficient, injection rate, and scattering time etc. 

Here, we use a dynamical Monte Carlo method to study the acceleration issue depending on factor 
of scattering time. In this isolated shock model, the maximum particle energy Emax will be calculated 
in different cases by applying different values for the constant of the scattering time. Since the FEB 
measures the size of the faded turbulent magnetic field in shock precursor region, if the FEB size is 
larger, then the Emax is higher. Unfortunately, if the size of the FEB is larger, then the computational 
expense is higher. Instead, we can change the scattering time to achieve a higher Emax in the shock. 
Assuming a particle can obtain the same additional energy gain from each cycle in a period of the 
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Table 1 Six cases with corresponding constants of the scattering time 


Simulation cases 






F 

t he scattering time 

ro 

to/2 

ro/'i 

ro/4 

to/5 

ro/12.5 


scattering time, it is probable that the more scattering probabilities will obtain more energy gains. If we 
take a smaller value for the constant of the scattering time in one simulation case, we can obtain a higher 
Emax value by more scattering probabilities in total simulation time. 

Although there still exists the impact from the diffusive coefficient in different shock region, it 
can be neglected in this isolated shock model. Since the ejecta moves with a large speed and the most 
part of accelerated particles are located in the turbulent magnetic field in the vicinity of ejecta, so the 
diffusive processes always can be taken as the Bohm condition and its coefficient difference in this 
limited precursor region would be slight. 

Monte Carlo approach regards the fluid as being composed of particles and focuses on the scatter¬ 
ing micro processes between the particles and turbulent magnetic field in the diffusion processes. This 
technique is based on the computational grids, where large number of particles distributed. Particle’s 
mean free path is proportional to its local velocity in its local frame as follows. 

\ = vl-t ( 1 ) 

where, t>L is the local velocity of part icles, r is the scat tering time. In Earth’s bow shock model, the 
scattering time r is taken as a constant dKnerr et ^Il996h . For comparing values for maximum particle 
energy Emax, we apply different values for constant of the scattering time to perform these correspond¬ 
ing cases as in Table [T] To simulate the scattering processes accurately, the scattering time r should be 
chosen to be far more than time step dt as follows. 

T :s> dt (2) 

To simulate the shock’s formation and evolution, we set the standard scattering time tq is a constant 
for all particles in Case A. Other constants of the scattering time and correspond ing cases can be see n 
in Tahiti] Those related simulation parameters can be refereed to previous work dWang & YanLl201^ . 
Here, we just list the scattering times in different cases. All of the scattering times are chosen to be more 
than time step dt (df=To/25) in corresponding cases. 

3 RESULTS 

3.1 Acceleration of particles 

To inspect the entire particles’ acceleration processes in the shock, we extracte a number of particles 
from the simulated box in each case. Six plots in Figure [T] are taken from six simulate cases from Cases 
A, B, C, D, E, to Case F, respectively. Each curve in each plot represents one particle’s evolution of 
velocity and energy with time. Every plot have a few peak velocities in some accelerated particles, 
the highest peak value in the velocity or energy axis represents the maximum velocity or energy in 
corresponding case. Each maximum energy value is denoted in each plot. Among of these six cases, the 
maximum energy value in Case C achieves an energy saturation at 5.5506MeV. In addition, we can also 
find that a part of particles at the bottom of each plot exhibit no acceleration in total simulation time. 
Another part of particles with jumps from the lower energy to the higher energy in each plot indicate 
their accelerating processes in the shock with time. Simultaneously, energy jumps in corresponding case 
show an increasing steeply tendency with a decreasing value for the constant of scattering time from 
Cases A, B, C, D, E, to Case F, respectively. Case F shows severely steep jumps and steeply descents 
in energy or velocity curves,because the scattering time r = ro/12.5 is chosen to approach to time step 
{dt = ro/25). These results indicate that the computational accuracy requires the scattering time should 
be more enough than time step. 
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Fig.l A number of particles are extracted from the simulation box in corresponding cases from Cases A,B,C,D, E, to Case F, respectively. Blue curve 
represents one particle's evolution of velocity and energy with time. The top peak in each plot shows the maximum velocity or energy in corresponding case. 
A part of particles have no acceleration at the bottom of each plot. A few curves with jumps from the lower energy to the higher energy in each plot indicate 
that they are accelerated in the shock with time. The Emax in Case C with a scattering time r = tq /3 achieves an energy saturation at 5.5506MeV 
with comparing six cases. 


3.2 The Emax Function. 

Here, we focus on an isolated CME-driven shock for calculating the maximum particle energy Emax 
in those cases applying the different values for the scattering time. Using our dynamical Monte Carlo 
model, we have obtained the different values for the Emax in those cases. So we can build the function 
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of the maximum particle energy Emax versus the values for the scattering time r with the value from 
To, ro/2, tq/3, tq/4, and tq/5 to to/12.5 in corresponding case. 


Utilizing the method described in Section Method, the calculations of Emax are performed un¬ 
der the scattering angular distribution with a standard deviation value for a = tt and an aver¬ 
age value for p = 0, which would relatively be mor e efficient for the p article acceleration in the 
CME-driven shock demonstrated by previous model (IWang & Yanl l2012h . From the large popula¬ 
tion of the accelerated particles at the end of the simulation in each case, we find each maximum 
local velocity VLmax in corresponding case with its scattering time r. The relationship between 
the maximum local velocity VLmax and the value for the scattering time r in all cases are shown 
in FiglJ] The solid line denotes the correlation of the maximum local velocity VLmax versus the 
different value for the scattering time r in corresponding case. In present Monte Carlo model, the 
value for the scattering time r is chosen from tq, tqI2, To/3, tq/T and tq/5 to to/12.5, respectively. 
And these six squares represent the maximum local velocity values for VLmax in all cases with 
their correspondingly values for the scattering time. The maximum local velocity VLmax is repre¬ 
sented by a dimensionless value for ULmax(yi)=17.7824, UL^a3.(s)=17.9928, UL^(j2,(c')=24.2773, 
VLmax{D)=i9-^'2-95, ULmax(£;)=17.3596, and VLmax(F)=^9.4-4S2 in each case, respectively. As 
shown in FigE] among these maximum local velocities VLmax, the largest is the one in Case C with 
a value of ULmo2;(c)=24.2773, and its value for the scattering time is tq/3. The top of the stair line in 
FigEJshows that there exists a saturation in a function of the maximum local velocity VLmax versus the 
value of the scattering time r under the resonant diffusion scenarios. 


Figl3]shows the fitting curve of the maximum particle energy Emax versus the values for the scat¬ 
tering time T. The maximum particle energy Emax are calculated in the shock frame by a scaled value 
according to the scale factor for velocity vscale - The maximum particle energy Emax in each case varies 
along the shape-preserving curve with a sequence of Emax(A)='^-92%GMe,V, £'max(B)=3.0489MeV, 
^^max(C)=5.5506MeV, £’„aa:(D)=3.7031MeV, F;„aa:(£)=2.8381MeV, and Emax{F) =3.5620MeV from 
Cases A, B, C, D, and F, to Case F. respectively. All of those maximum particle energy Emax are not 
exceed to the upper limit of F^reafe at lOMeV in the observation. But Case C with corresponding value 
for the scattering time tq/S shows that the largest maximum particle energy is £'maa:(C)=5.5506MeV, 
which is still less than the upper limit of the E^reak region. It implies that whatever the value for the 
scattering time is chosen under an isolated shock model, the maximum particle energy Emax would not 
more than the upper limit of the Ebreak region in the observed energy spectrum. According to these 
simulation results, the energy spectrum “cut-off’ would be formed near the energy range of 5MeV. 
Simultaneously, the saturation value for the maximum energy function demonstrate that these maxi¬ 
mum particle energy Emax can fit the observed lower energy spectrum below the E^reak limit. Booking 
from the shape-preserving curve in Figl3 the Emax will not increase with the value for the scattering 
time T decreasing from tqI3 to tq/S. Although the function of maximum particle energy Emax shows a 
lightly ascendant tendency when the value for the scattering time decreasing from tq/S to ro/12.5, the 
value for the scattering time ro/12.5 is approaching to the time step dt (i.e., dt = to/25). Considering 
the precision of the calculation, the value for the scattering time t should be chosen to be not less than 
the time step dt. Since the amplified magnetic field is limited by the order of the magnitude SB/B ^ 1, 
whatever the value for the scattering time is chosen in an isolated shock model, the obtained maximum 
particle energy Emax are not more than the upper limit of the Ebreak in the observed energy spectrum. If 
we expected to obtain a more extended energy spectrum beyond the upper limit of Ebreak at lOMeV and 
even up to GeV, the multiple shocks model would be applied. So it means that the efficiency of acceler¬ 
ation in an isolated shock model will not exceed the upper limit of the Ebreak as long as the value for 
the scattering time is chosen to be enough more than the time step dt. In addition, it also implies that the 
realistic observation of the Ebreak energy spectrum strongly requires a multiple shock model to transfer 
the shock’s energy into the superthermal particles up to a highest energy spectrum for explaining the 
Ebreak formation and the higher energy spectrum. 






ENERGY SPECTRAL PROPERTY IN A CME-DRIVEN SHOCK 


7 


30 

28 

(n 

o 

I 26 
(0 

I 24 

'cj 

^22 
> 

w 20 
o ” 
o 

_i 

P 18 


12 

10 


Comparison of the VLmax in Different Cases 


Cases A 


B 


24.2773 


19.8295 


19.4482' 


17.3596 


17.7824 17.9928 

I ■ Simualtion Cases in A,B,C,D,E,F(g=7t,^=0) | 


1/2 1/3 1/4 

The Scattering Time T (x ) 


1/5 2/25 


Fig. 2 The maximum local velocities VLmax of the accelerated energetic particles in all cases 
are plotted as a stair function of the value for the scattering time. All these cases are simulated 
using Monte Carlo model complying with a Gaussian scattering angular distribution with a 
standard deviation value for a = tt, and average value for /i = 0. 


Fitting Function of the Maximum Energy in Cases with Different t 



Fig. 3 The shape-preserving curve of the maximum particle energy Emax as a function of 
the value for the scattering time in six cases. The values for maximum particle energy are 
indicated by scaled values along the profile of the shape-preserving curve in all cases with 
different value for the scattering time. The profile of the Emax in different cases shows a 
saturation value for the maximum particle energy in Case C with a critical value for the 
scattering time to/3. But all of them are not exceed the upper limit of the E^reak at 10 MeV 
of observation. 


3.3 The Energy Spectra 

FiglDshows the shock energy spectra calculated in the downstream region in all cases. As far as the shape 
of the energy spectrum is concerned, the power-law slope of six extended curves are similar, because all 
cases are done in the same resonant diffusion scenarios just only with different values for the scattering 
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time. However, among these cases, the energy spectrum in Case C with the value for the scattering time 
of to/ 3 shows a relatively hard slope in the highest energy spectral tail. Under an isolated shock model, 
each case shows that how the initial Maxwellian energy spectrum to evolve into the extended energy 
spectrum with “power-law” structure in its high energy, respectively. By comparison, we calculated 
the average value of the maximum particle energy in present six cases. The average value for maximum 
particle energy is < Emax >=3.6135MeV and the average value for energy spectral index is F ~ 1.125. 
These results agree with t he low energy spectr um in the observations from the multiple spacecraft. 
Observed energy spectrum dMewaldt et al.Ll200^ shows low energy spectrum with an index of F = 1.07 
and a high energy spectrum with an index of F = 2.45. The observed energy spectrum indicated that 
there exists an Ebreak between the lower energy spectrum and the higher energy spectrum. From these 
simulated cases, we concluded that all these energy spectra are characterized by a “power-law” with 
an averaged index F ^ 1.125, which consists with the observed index F = 1.07 of the low energy 
spectrum. Since there is no maximum particle energy Emax in these six cases beyond the upper limit of 
Ebreak at lOMeV, we are unable to conclude that there should exist an Ebreak at 1-lOMeV as a “break” 
from the lower energy spectrum to higher energy spectrum at this range. If we expect the second higher 
energy spectrum to occur, we can guess that there must exist an enhancement in amplification of the 
magnetic field associated with multiple shocks model. We propose the multiple shocks model would 
be applied to further investigate the higher energy spectrum in CME shock events. Recen tly, there are 
some analysis of multiple CME collision events that have been discussed. Eor example, ICheng et al.l 
(l2013h report the initiation process of compound CME activity consisting of t wo succes s ive er uptions 
of flare ropes that occurred on January 23 2012. Another example presented bv iLiu et"^ (l2014l) shows 
that the interactions between consecutive CMEs resulted in a “perfect storm” near 1 AU on 23 July 
2012, which would induce to nonlinear amplification of magnetic field. Eurther more evidences could 
be gathered from observations in spacecraft such as SDO, SOHO, ACE, Wind and etc. In the implication 
from these present simulated results, we propose to build a multiple shocks model to simulate the Ebreak 
formation and the higher energy spectrum in the interplanetary shock. In present model, we think that the 
parameter of the scattering time would play key role on the strength of the diffusive coefficient for Emax 
production within the resonant diffusion scenarios associated with isolated shock. According to the final 
results, we find the relationship between the maximum particle energy Emax and the different value for 
the scattering time in isolated shock model. Although the difference on these maximum particle energy 
Emax in simulated cases has happened, no maximum particle energy Emax can exceed the upper limit 
of Ebreak to further evolve into a higher energy spectrum up to GeV. The future simulation is necessary 
to verify the higher energy spectrum with an index of F ~ 2.5 and the energy spectral “break” formation 
by applying multiple shocks model. 

4 SUMMARY AND CONCLUSIONS 

In summary, these presented simulations are unable to verify that there should exist an energy spectral 
“break” below lOMeV in some large CME-driven shocks. But instead we verify that there is an energy 
spectral “cut-off’ near the range of energy at 5MeV in an isolated CME-driven shock. We calculate 
the maximum particle energy Emax focusing on Dec 14 2006 CME-driven shock event and built the 
relationship between the maximum particle energy Emax and the value for the scattering time r. We find 
that the maximum particle energy Emax approaches to a saturation near 5MeV below the upper limit 
of Ebreak of the observed energy spectrum. We verify the lower energy spectrum is consistent with the 
observed low energy spectrum, but no higher energy spectrum appears. Although there are several large 
SEP events in the past solar cycle 23 appear the energy spectral “breaks” between 1-lOMeV, there is 
still no very reasonable explanation. Since these observations are depend on multiple spacecraft, it is not 
easy to treat the system errors and couple the observed data obtained from different spatial orientations. 
The huge computational expense also limits numerical methods to reach an enough high energy spectral 
tail for further to identify this “break”. In the view of the current theoretical point about DSA, analytic 
method give an implication that this “break” would be connected with particle leakage mechanism. 
This “break” seemly can be predicted in the escaped position ahead of supernova remnants (SNRs) 
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Fig. 4 The energy spectra obtained from downstream region in six cases with different values 
for the scattering time. The thick solid line with a narrow peak at E =1.4315keV represents 
the initial Maxwellian energy distribution in the upstream region. All of cases consist with 
the low energy spectrum with all Emax less than lOMeV. 


shock, where the SNRs has a collision with nearby molecule clouds (MC). This idea will motivate us 
to further investigate the energy spectrum E^reak formation and the higher energy spectrum. Hopefully, 
we propose the multiple shocks model would play key role on explaining the energy spectrum E^reak 
formation and the higher energy spectral shape. 
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